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ABSTRACT 

We study various measures of weak lensing distortions in future surveys, taking into 
account the noise arising from the finite survey size and the intrinsic ellipticity of 
galaxies. We also consider a realistic redshift distribution of the sources, as expected 
for the SNAP mission. We focus on the low order moments and the full distribution 
function (pdf) of the aperture- mass M ap and of the smoothed shear component "fi S . We 
also propose new unbiased estimators for low-order cumulants which have less scatter 
than the usual estimators of non-Gaussianity based on the moments themselves. Then, 
using an analytical model which has already been seen to provide a good description 
of weak gravitational lensing through comparison against numerical simulations, we 
study the statistical measures which can be extracted from future surveys like the 
SNAP experiment. We recover the fact that at small angular scales (1' < 9 S < 10') 
the variance can be extracted with a few percent level accuracy. Non-Gaussianity can 
also be measured from the skewness of the aperture-mass (at a 10% level) while the 
shear kurtosis is more noisy and cannot be easily measured beyond 6'. On the other 
hand, we find that the pdf of the estimator associated with the aperture-mass can 
be distinguished both from the Gaussian and the Edgeworth expansion and could 
provide useful constraints, while this appears to be difficult to realize with the shear 
component. Finally, we investigate various survey strategies and the possibility to 
perform a redshift binning of the sample. 

Key words: Cosmology: theory - gravitational lensing - large-scale structure of 
Universe - Methods: analytical, statistical, numerical 



1 INTRODUCTION 

Weak lensing surveys have already started playing a major 
role in cosmology (e.g., Bacon, Refregier & Ellis, 2000, Hoek- 
stra et al., 2002, Van Waerbeke et al., 2001, and Van Waer- 
beke et al., 2002) not only in constraining the background 
dynamics of the universe but in probing the nature of dark 
matter as well. In order to extract useful information from 
observations one needs to compare the data with theoretical 
predictions associated with specific cosmological scenarios. 
To do so, one often uses numerical simulations which typi- 
cally employ ray-tracing techniques as well as line of sight 
integration of cosmic shear (e.g., Schneider & Weiss, 1988, 
Jaroszynski et al., 1990, Wambsganss, Cen & Ostriker, 1998, 
Van Waerbeke, Bernardeau & Mellier, 1999, and Jain, Seljak 
& White, 2000, Couchman, Barber & Thomas (1999)). On 
the other hand, several analytical techniques have also been 



developed over the past several years to predict statistics of 
weak lensing shear and associated quantities. On large angu- 
lar scales perturbative techniques are generally employed as 
non-linearities can be treated by a series expansion (e.g., 
Villumsen, 1996, Bernardeau et al., 1997, Jain & Seljak, 
1997, Kaiser, 1998, Van Waerbeke, Bernardeau & Mellier, 
1999, and Schneider et al., 1998). However, on small angular 
scales, especially relevant to current observational surveys 
with small sky coverage, perturbative calculations are no 
longer valid and models to represent the gravitational clus- 
tering in the non-linear regime have had to be devised. One 
class of such models is based on the hierarchical ansatz (e.g., 
Fry 1984, Schaeffer 1984, Bernardeau & Schaeffer 1992, and 
Szapudi & Szalay 1993, 1997, Munshi, Melott & Coles 1999) 
for the evolution of high-order correlation functions, joined 
with Peacock & Dodds (1996)'s prescription (see Peacock & 
Smith (2000) for more recent fit) to model the evolution of 
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the power spectrum or equivalently the two-point correla- 
tion function. Applications of such hierarchical models have 
shown to be quite accurate in predicting various statistics 
related to weak lensing shear and convergence at small an- 
gular scales (Valageas 2000a & b; Munshi & Jain 2000 & 
2001; Munshi 2000; Bernardeau & Valageas 2000; Valageas, 
Barber & Munshi 2004; Barber, Munshi & Valageas 2004; 
Munshi, Valageas & Barber 2004). A second class of theo- 
retical descriptions of the density field is based on the halo 
models (see Cooray & Seth 2002 for a review) which can also 
reproduce lower order moments (e.g. Takada & Jain 2002, 
Takada & Jain 2003a,b). 

In order to make contact with observations, in addition 
to a good description of the underlying matter density field, 
which gives rise to these weak lensing effects, one clearly 
needs to include various sources of noise such as the contri- 
butions from the intrinsic ellipticity distribution of galaxies, 
shot-noise due to the discreet nature of the source galaxies 
and finite volume effects due to finite survey size. Following 
Schneider et al. (1998) a detailed analysis of these effects 
was presented in Munshi & Coles (2003). We extend these 
results in some respects by incorporating a realistic redshift 
distribution of sources and generalizing the estimators used 
to directly handle the components of shear. We also propose 
a new family of estimators for low order moments which have 
less scatter than usual ones. We study in details the influ- 
ence of the source redshift distribution and of the intrinsic 
galaxy ellipticities on the measures, as well as the sensitiv- 
ity to various cosmological or survey parameters. Extending 
our analysis to complete probability distribution functions 
we investigate to what extent various noise contributions 
make it difficult to distinguish the signatures of the under- 
lying non-linear dynamics hidden in the tails of the pdf or 
its shape near its maximum. In most cases we use for nu- 
merical display the parameters associated with the SNAP 
experiment. 

This paper is organized as follows: in section 2, we de- 
scribe weak lensing observables in general and various filters 
used to smooth the data, as well as the estimation of the ob- 
servational scatter. In section 3, we introduce specific survey 
geometries based on SNAP class experiments and we com- 
pute the noisy cumulants and the associated probability dis- 
tribution functions for estimators of the shear components 
as well as the aperture mass. Finally, in section 4 we discuss 
our results. 



2 STATISTICS OF WEAK-LENSING 
OBSERVABLE 

In a series of papers (Valageas et al. 2004, Barber et al. 2004, 
Munshi et al. 2004) we have described how to obtain the dis- 
tribution function of any weak-lensing observables and we 
have shown that our predictions match results from numer- 
ical simulations for the specific cases of the convergence k, 
the shear 7 and the aperture- mass M ap . However, in those 
studies we did not include the noise associated with the in- 
trinsic ellipticity of galaxies and we assumed that all sources 
were located at the same redshift z s - Here, we show that our 
formalism can be extended in a straightforward manner in 
order to handle these two effects. 



2.1 Redshift distribution of sources 

Let us first recall our notations. Weak-lensing effects can be 
expressed in terms of the convergence along the line-of-sight 
towards the direction i? on the sky up to the redshift z s of 
the source, k(i9,z s ), given by (e.g., Bernardeau et al. 1997; 
Kaiser 1998): 



*,{•&, z s ) = 



with: 



3f2 n 
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(1) 
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where z corresponds to the radial distance x an d T> is the 
angular distance. Here and in the following we use the Born 
approximation which is well-suited to weak-lensing studies: 
the fluctuations of the gravitational potential are computed 
along the unperturbed trajectory of the photon (Kaiser 
1992). Thus the convergence k(-&, z s ) is merely the projec- 
tion of the local density contrast S along the line-of-sight. 
Therefore, weak lensing observations allow us to measure the 
projected density field k on the sky (note that by looking 
at sources located at different redshifts one may also probe 
the radial direction). In practice the sources have a broad 
redshift distribution which needs to be taken into account. 
Thus, the quantity of interest is actually: 



k(#) 



f 

Jo 



dz s n(z s )ft(i9, z 3 ), with / dz s n(z s ) = 1, (3) 



where n(z s ) is the mean redshift distribution of the sources 
(e.g. galaxies) normalized to unity. From eq.(l), the conver- 
gence k associated with a specific survey also reads: 



k(i?) = 



with: 
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dx w(x) <5(x,2>#), 



dz 3 n(z s ) w(\,Xs), 



(4) 



(5) 



where z max is the depth of the survey (i.e. n(z s ) — for 
z s > 2 max ). By working with eq.(4) we neglect the discrete 
effects due to the finite number of galaxies. They can be 
obtained by taking into account the discrete nature of the 
distribution n(z s ). This gives corrections of order 1/N to 
higher-order moments of weak-lensing observables, where N 
is the number of galaxies within the circular field of interest. 
In practice N is much larger than unity (for a circular win- 
dow of radius 1 arcmin we expect N J> 100 for the SNAP 
mission) therefore in this paper we shall work with eq.(4). 

Thus, in order to take into account the redshift distribu- 
tion of sources we simply need to replace w(x, Xs) in eq.(l) 
by w(x)- Therefore, all the results of Munshi et al. (2004) 
remain valid. Then, usual weak-lensing observables can be 
written as the angular average of /t($) with some filter U: 



I 



X = / dtf Ux(0) k(i9). 



(6) 



For instance, the filters associated with the smoothed con- 
vergence k s , the smoothed shear 7 S and the aperture- mass 



M a p are (Munshi et al. 2004): 



9(fl < 9s) 
n9 2 s ' 



e(fl > 6 S ) l2a 

TTt? 2 6 



(7) 
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and 



(8) 



where O are Heaviside functions with obvious notations and 
q is the polar angle of the vector The angular radius 
8s gives the angular scale probed by these smoothed ob- 
servables. Note that the smoothed shear 7 S depends on the 
matter located outside of the cone of radius 8 S . However, in 
practice one directly measures the shear 7(1?) on the direc- 
tion ■& (from the ellipticity of a galaxy) and 7 S is simply the 
mean shear within the radius 8 S . For M ap we shall use in this 
paper the filter (8), as in Schneider (1996), but one could 
also use any compensated filter with radial symmetry. As 
in Valageas (2000a), it is convenient to define the minimum 
convergence K m j n associated with an empty beam (5 = — 1): 



" 2 J 



and to normalize all observables with respect to n n 
X 



X = 



with: 



d X w / MUx(&) 5{x,m), 



(9) 



(10) 



(11) 



Then, as described in Munshi et al. (2004), the cumulants 
of X can be written as: 



X 



dx w p 




X2 

vS 2 



Xp 



(12) 



J i=l 



dk ±j Wx(k ±j V8 s ) 



x (<y(k±i)...J(k±p)> t 



(13) 



We note (...) the average over different realizations of the 
density field, £ p is the real-space p— point correlation func- 
tion of the density field £ p (xi, . . . , x p ) = (<S(xi) . . . <5(x p )) c , 
fc|| is the component of k parallel to the line-of-sight and kx 
is the two-dimensional vector formed by the components of 
k perpendicular to the line-of-sight. In eq.(13) we factor- 
ized the Dirac term 8o{kui + . . . + k^ p ) out of the connected 
correlation (5(k±i) . . . 8(k± p )) c . We note Wx(k±T>8 e ) the 
Fourier transform of the window Ux: 



Wx(k±V8 3 ) = / dtiUxW) 



ik 1 .X>t? 



(14) 



In particular, for the smoothed convergence k s , the smoothed 
shear 7 S and the aperture-mass M ap we have (Munshi et al. 
2004): 



W Ks (k ± V8 s ) = 



2J 1 (k ± V8 s ) 
k ± V8 s : 



and using eq.(8): 



(15) 
(16) 



W M ^(k±V8 s ) 



24J 4 {k ± V8 s ) 
(k ± V8 s y ■ 



(17) 



where <f> is the polar angle of kx and J v are Bessel functions 
of the first kind. The real-space expression (12) is well-suited 
to models which give an analytic expression for the cor- 
relations £ p , like the minimal tree-model (Valageas 2000b; 
Bernardeau & Valageas 2000; Barber et al. 2004) while the 
Fourier-space expression (13) is convenient for models which 
give a simple expression for the correlations (<5(ki)..<5(k p )) c , 
like the stellar model (Valageas et al. 2004; Barber et al. 
2004). In these two cases one can resum the cumulants {X p ) c 
which yields the pdiV(X) as (e.g., Munshi et al. 2004): 



V(X) 



d V \Xy- Vji (y)}/(X 2 ) 



■ v. 2ni{X*} c 
where we introduced: 



(18) 



(19) 



The generating function tPx{y) is closely related to the char- 
acteristic function <p(y) of the density field. Thus, for the 
stellar model <p x (y) is merely a suitable average along the 
line-of-sight of ip{y), while for the minimal tree-model the re- 
lationship is slightly more intricate but explicitly known. For 
the smoothed convergence, one actually has iph e {y) — f(y) 
(Valageas 2000a,b; Barber et al. 2004). 



2.2 Intrinsic ellipticity of galaxies, pdf of 
weak-lensing observables 

2.2.1 Aperture-mass M ap 

Th expressions obtained in the previous section assumed 
that the observations were perfect. However, in practice the 
data exhibits some noise. A specific source of noise is merely 
due to the intrinsic ellipticity of galaxies, which cannot be 
avoided. Thus, in order to measure the aperture-mass M ap 
within a single circular field of angular radius 8 S , in which 
iV galaxies are observed at positions i9j with tangential el- 
lipticity e t ,j, we can use the estimator M defined by: 



M 



:*f ap j e%,j- 



(20) 



j=i 



Here we used the fact that the aperture-mass defined from 
the convergence k by the compensated filter Um^ given in 
cq.(8) can also be written as a function of the tangential 
shear 7t as (Kaiser et al. 1994; Schneider 1996): 



M a 



with: 



/ 



dtf Q M ap ($) 7tC#); 



Qm, p (0) 



< 



f- 

\8s 



fl2 



(21) 



(22) 



In eq.(20) we wrote Qj\ 



Qm 1p In the case of weak 



lensing, k< 1, the observed complex ellipticity e is related 
to the shear 7 by: e — 7 + e» , where e* is the intrinsic ellip- 
ticity of the galaxy. Assuming that the intrinsic ellipticities 
of different galaxies are uncorrelated random Gaussian vari- 
ables, the cumulant of order p of M is: 
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(M p ) c = <Mf p ) c 1 + 



2iV(M a 2 p ) 
Wlth P= 2 n ' ( 23 ) 

where (5 Pj 2 is the Kronecker symbol, a» = (e») is the disper- 
sion of the intrinsic ellipticity of galaxies and we introduced: 



7 



/ dtfQ Map W 2 



(24) 



For the filter (22) we have G?M ap = 6/5. In order to ob- 
tain eq.(23) we have averaged i) over the intrinsic ellipticity 
distribution, ii) over the galaxy positions and iii) over the 
matter density field, assuming these three averaging proce- 
dures are uncorrelated. The second step can be written for 
any quantity X as: 



N 

(^)galaxy positions = J^J / dZj Tl^Zj^ I 



Mi y 



(25) 



a(P k f = (Pi) - (P k ) 2 = g 



P k A 



(31) 



Of course, we recover the scaling a(P k ) oc 1/-/N^ where N c 
is the number of cells. On the other hand, since different 
intervals I k do not overlap their covariance is: 



fc/fc': a(P k ,P k ,) 2 = (P k Py) 



PkPy — -j-j-PkPk 1 - (32) 



2.2.2 Smoothed shear component f iB 

In a similar fashion to the aperture-mass, we can measure 
the shear-component 7* (with i = 1,2) from the estimator 
Ti which we define by: 



- N 2^ 



(33) 



Thus, since the intrinsic ellipticities are Gaussian and we 
neglected any cross-correlation with the density field they 
only contribute to the variance of the estimator M. Note 
that M 2 is a biased estimator of (M 2 p ) because of this ad- 
ditional term. The quantity p measures the relative impor- 
tance of the galaxy intrinsic ellipticities in the signal. They 
can be neglected if p 3> 1. Note that any Gaussian white 
noise associated with the detector can be incorporated into 
the expression (23) by adding a relevant correction to a 2 . Fi- 
nally, from eq.(23) and eq.(19) we obtain for the generating 
function ip^ of the normalized quantity M — M/|/t m i n |: 



1+P 



1 + p 



i v 2 

l + p 2 



(26) 



Of course, for small p we recover the Gaussian (i.e. (j/) = 
— y /2) while for large p we recover ip M ~ (y). 

Thus, each circular field of angular radius 8 S yields a 
particular value for the quantity M defined in eq.(20). If the 
survey contains N c such cells on the sky, we can estimate 
the pdf V(M) through the estimator: 



P k = 



N C A ^ 

n = l 



Ufa), 



(27) 



where l k (n) is the characteristic function of the interval I k 
of width A, applied to the value M(n) of M measured in 
the cell n: 



lk{n) = 1 if M(n) £ I k , l k (n) = otherwise, 



with: 
h = 



A A 

M k --,M k + - 



M k = kA. 



(28) 



(29) 



For simplicity, we chose all intervals I k to have the same 
width A, but this could easily be modified. Note that dif- 
ferent intervals do not overlap and the integer index k runs 
from — oo to +oo. Then, from the set {P k } we obtain an his- 
togram which provides an approximation to V{M). Indeed, 
we have: 

-M k +A/2 



(Pk) = Pk With P k = 



dM 



V(M). 



(30) 



M fc -A/2 



Then, for small enough A we have Pk ^ P(M k ). Next, 
assuming that different cells are uncorrelated, the dispersion 
of the estimator P k is: 



Here eij is the component i of the ellipticity of the galaxy j 
and for the smoothed shear component 7i S we have Q 7i ($) = 
1/(tt9 2 ) which is independent of &. Thus we now get G lis = 1 
and we recover eq.(26) relating v? r -. (y) to <fi-y~ its (y), where we 
now use G lig = 1 for p (not to introduce too many notations, 
we use the same letter p for both the aperture-mass and the 
shear) . 

Next, we can estimate the pdf V(Ti) as in eq.(27). How- 
ever, we can take advantage of the fact that the pdf P(Ti) 
is even. Therefore, we can group the intervals and I k to 
evaluate P{T itk ). In other words, we now write: 



Pk 



and: 



iV c 2A 



with k — 0, 1, 2, .., oo, 



h 
where: 



fi.-fc - — ,r 



A 

- k + - 



li fe , Lik ~\~ 



k + 



This yields: 

(Pk) = Pk 



and 



p 2 

/ o \2 r k 

a{Pk) =w c 



P k 2A 



(34) 



(35) 



(36) 



(37) 



where P k is defined as in eq.(30). Thus, since V(Ti) is even 
we have gained a factor 2 in the expression (37) of the disper- 
sion a(P k ) 2 . On the other hand, for k / k' , the covariance 
is again given by eq.(32). 

Finally, let us note that we kept the term associated 
with the galaxy intrinsic ellipticities, which scales as 1/iV, 
while we neglected the terms associated with the fluctu- 
ations of the redshift and angular distribution of sources, 
which also scale as 1/N. The reason for doing so is that the 
correction due to the galaxy intrinsic ellipticities involves 
the multiplicative factor o~ 2 /{M 2 p } which can be large so 
that 1/p can be large even though we have N ^> 1. 



2.3 Low-order moments 

2.3.1 Aperture-mass M ap 

The quantities M and 7» introduced in the previous section 
provide biased estimators for the moments of weak-lensing 
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observables. In practice it is desirable to build unbiased esti- 
mators in order to measure low-order moments. Thus, as in 
Schneider et al. (1998) or Munshi & Coles (2003), in order 
to study the aperture-mass we can define the estimators M p 
as: 



Mr, 



2\p 



N(N — l)..(iV — p + 1) 



I) 51 < 3ji e tji---0jp e t, Jp , (38) 



The sum runs over all combinations {ji, .., j p } with no iden- 
tical indices. For simplicity, we wrote Q for Qm bp • If we are 
interested in the smoothed shear component 7; s we simply 
need to use Q 7i ,j and a t j in eq.(38). Then, a straightforward 
calculation gives the expectation values of the estimators M p 
as well as their dispersion a(M p ) 2 : 



(M p ) = (Ml p ), a(M p ) 2 = (M 2 ) - (M p y 
with: 

c(Mi) 2 = (M! p ) (l + -Y 



a(M 2 ) 2 = (M a 4 p ) c + 2{Ml p ) 2 1 + 



(39) 

(40) 
(41) 



a(M 3 ) 2 = <M a b p ) c + <M a 4 p ) c <M 2 p > ^15 + - j + 9<M aP .> 2 

+ «) 3 (l5 + ^ + i| + ^), (42) 

where p, which was defined in eq.(23), measures the contri- 
bution of the "cosmic variance" to the noise, relative to the 
galaxy intrinsic ellipticities (and detector white noise). Of 
course, the dispersion of M p involves the cumulants of M ap 
up to order 2p. In eqs.(40)-(42) we assumed N 3> 1 and we 
neglected relative corrections of order 1/N. The estimators 
M p correspond to a single circular field of angular radius 9 S 
which contains N galaxies. In practice, the size of the survey 
is much larger than 9 a and we can average over N c cells on 
the sky. Thus, we define the estimators M p as: 



Mi 



-T 



M v 



(43) 



where M p , n is the estimator M. 



for the cell n. Assuming 
that these cells are sufficiently well separated so as to be 
statistically independent, we have: 



(M p ) = (M p ) = (M*,), a(M P ) 



<t(M p ) 



(44) 



Here we assumed for simplicity that all cells have the same 
number N of galaxies. The skewness of the aperture-mass 
M ap is the coefficient S^ p defined as in eq.(19): S^ Up = 
<M a 3 p ) c /(M a 2 p ) 2 . Therefore, it can be estimated from the ra- 
tio 53: 



S 3 = 



Ms 



ajMs) 
(Ml) 2 ' 



(45) 



In eq.(45) we neglected the dispersion of M 2 in order to 
obtain the mean and the dispersion of 53. This is not a 
serious shortcoming because the error bars increase very fast 
with the order of the moments so that the dispersion of 
53 is dominated by the dispersion of M 3 . To avoid this 



approximation one may simply study the cumulants (M£ p ) c 
themselves rather than the ratios SiY^ . 



2.3.2 Smoothed shear component ji s 

For the smoothed shear component 7» s we obtain in a similar 
fashion: 



(M P ) = (M p ) = < 7 f s ), a(M P ) 
with: 



o-(M 2 ) 2 = < 7 f s ) c + 2< 7 2 ) 2 [l + - p 



a(M p ) 



(46) 



(47) 



a(M 4 ) 2 



(-:>. + :-!:^.(- 2 : I 28+^ ) +:iJ(- 



+(7i>c(7f s ) 2 204 + 



+dlr 96 + 



240 72 



240 216 96 24 



P 



2 + — + 

pj. pA pi 



(48) 



Since the shear component 7, is an even quantity (its sign 
can be changed by a simple change of coordinates) all odd 
moments vanish. The kurtosis of the shear SJ" is defined as 
in eq.(19): 



SJ' 



<7^)c _ <7f s > " 3(7£> i 



<7 2 s ) s <7 2 s ) s ' 

Therefore, it may be estimated from the ratio 54: 



5 4 = 



M% 



*J4 1 



a(5 4 ) 



(49) 



(50) 



<7 2 > 3 

Here we again neglected the dispersion of M 2 - We must 
point out that while the correlation between different cells 
on the sky becomes quickly negligible as soon as they do 
not overlap if we consider the aperture-mass (Schneider et 
al. 1998) this is not the case for the shear itself. This is due to 
the fact that M ap involves a compensated filter which damps 
the contribution from long wavelengths {Wm^ v (k±T>9 s ) ~ 
(k ± V9 s ) 2 for k ± V6 s < 1 in eq.(17)), while the shear is 
sensitive to low k (|W / 7s (k_ L r'6' s )| -> 1 for k -> in cq.(16)). 
In this article we shall not investigate this point and we shall 
assume that the N c cells are sufficiently far apart so as to 
exhibit a negligible correlation. 

2.4 Improving low-order estimators 

The interest of the cumulant (Mf p ) c is that it provides a 
measure of the deviations from Gaussianity. Moreover, it 
can be used to break the degeneracy between the cosmolog- 
ical parameter O m and the normalization ag of the matter 
power-spectrum (Bernardeau et al. 1997). However, we can 
note that the estimator M3 may not be the best way to mea- 
sure the skewness. Indeed, in order to make the most of the 
departure of the pdf V(M ap ) from the Gaussian, we would 
like to weight the pdf by a factor which changes sign with 
the difference - P(M ap ) — 'pG(Afap), where Vg is the Gaussian 
with the same variance as "P. A detailed study of V(M Bp ) 
was presented in Munshi et al. (2004), see also section 3.1.3 
below. It shows that ^(Map) — Pa(M ap ) does not change 
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sign with M ap (like Mf p ). On the other hand, for small de- 
viations from Gaussianity the pdf of a quantity X can be 
expanded as (e.g., Bernardeau & Kofman 1995): 



V{X) = 



1 



- v 2 / 2 



24 



72 



H e {u) 



with: 



ax = (X 2 ) 1 ^ 2 and v = 



X_ 

ax 



(51) 



(52) 



Here we introduced the Hermite polynomials H n (y). In par- 
ticular we have: 



H-.i(v) — v 3 — 3v and Hi(y) — u* — &u 2 + 3. 



(53) 



The Edgeworth expansion (51) is only useful for moderate 
deviations from the Gaussian, that is when the first cor- 
recting term is smaller than unity (typically \u\ Si 1 and 
\(?xSz\ £ 1). As seen in Munshi et al. (2004), the Edge- 
worth expansion (51) is actually not very useful to describe 
the pdf of the aperture-mass or the shear, since it only fares 
well when the pdf is very close to Gaussian. However, as we 
have seen in section 2.2, the noise due to the intrinsic ellip- 
ticity of galaxies makes the pdf of actual observables closer 
to Gaussian. Moreover, even if the pdf obtained from eq.(51) 
is not very accurate it gives a reasonable prediction for the 
sign of the difference V(M a , p ) — Vg(M &p ). This leads us to 
consider for the case of the aperture-mass the estimators Hz 
and Hz defined by: 



H 3 = M3-3M 2 M!, 



1 

Hz = — H 3 , n , 



(54) 



where M p and M v are the estimators introduced in sec- 
tion 2.3. Thus, Hz involves the estimators M3 and Mi as- 
sociated with the same cell, as well as the mean M 2 over 
all cells. Of course, H3 is built from the Hermite polynomial 
Hz(v) given in eq.(53). We see from eq.(54) that it is actu- 
ally an estimator for the cumulant (M ap ) c rather than for 
the moment (M 3 p ) like M3 (but in the case of the aperture- 
mass it happens that (M 3 p ) c = (M 3 p )). At lowest order over 
1/N C we obtain: 



(Hz) = (Hz) = «,)„, a(Hz) ~ a{Ha) 



with: 



a(H 3 f = *(Mz) 2 -6<M 4 p } c <M a 2 p > -9<M a 2 p > 3 1 + 



(55) 



.(56) 



Therefore, we see that the dispersion of Hz is always smaller 
than for A4z- Note that although Hz as written in eq.(54) 
is biased by a term of order 1/N C , which should not be a 
problem, this term can be removed in a straightforward way 
by replacing M 2 in eq.(54) by (N C M 2 - M 2 , n )/(N C - 1) (i.e. 
the mean M 2 is computed from all other (N c — 1) cells). 
Then, the skewness S^ 1 ' can be estimated from 

<t(Hz) 



Q H _ Hz 
°3 ~ T72 



<5 ; f) 



(m 2 p ) 2 ' 

In eq.(57) we again neglected the dispersion of M 2 . 



(57) 



For the shear, we introduce in a similar fashion the es- 
timators H4 and H4 built from the Hermite polynomial of 
order four: 



H 4 = M 4 - 6M 2 M 2 + 3Mi, Ha 



-y 



Hi 



(58) 



This again provides an estimator for the cumulant (-yt B ) c 
rather than the moment = (7, 4 S ) C + 3(7, 2 S ) 2 which was 
estimated by M4. Then we obtain at lowest order over 1/N C : 



(H4) = (H 4 ) 
with: 



(lt)c, <y(Hi) = 



°(Hi. 



(59) 



a(Hi) 2 = a(Mi) 2 - 12( 7 f s ) c ( 7 2 ) 



n + 8 - 

p 



72< 7 ? s ) 4 



l + i 

P 



(60) 



Thus, we note that the dispersion of Hi is again always 
smaller than for Mi- Next, the kurtosis SJ" may be esti- 
mated from: 



Sf = 



Hi 



<5f> 



a {Hi) 

(ll) 2 ' 



(61) 



where we again neglected the dispersion of M 2 . 

Thus, since the estimators Hz and Hi are no more dif- 
ficult to use than M3 and M4 and exhibit a smaller dis- 
persion, they are better suited to the measure of low-order 
cumulants. In fact, it can be shown that Hz and Hi are 
optimal estimators among their class for a Gaussian distri- 
bution. Thus, if we generalize Hz and Hi to the unbiased 
estimators Lz and Li defined as: 

L 3 = M 3 -a 3 M 2 M 1 , L 4 = Mi-aiM 2 M 2 +(ai-3)Ml,(62) 

where az and 04 are free parameters, one can easily see that 
(Lz) = (M 3 p ) c and (Li) = (7j 4 s ) c and that the variance of 
these estimators is minimum for : 



Q3 = 3 + 



and: 



<M a 4 p ) c 



<Mf p ) 2 (l + l/p)' 



(63) 



a.i 



_ 12( 7 2 ) 4 (1 + 1/p) 2 + (li) c (ll) 2 (U + 8/p) + (7|)cv7|) 

o/-,2\4/1 1 1 I „\ 2 1 /».4\ /„,2\2 



2<7 2 ) 4 (1 + 1/P) 2 + <7 4 ) C (7 2 ) 2 



(64) 



Therefore, for a Gaussian distribution the estimators Lz and 
Li show the lowest scatter for az — 3 and a 4 = 6, that 
is when they are identical to the estimators Hz and Hi. 
Since weak lensing observables are not exactly Gaussian, 
using L p instead of H v could further lower the dispersion of 
the measures. However, since the cumulants are not known 
a priori (these are the quantities to be measured !) it is 
probably better to use the simple estimators H v . Moreover, 
as seen above the estimators H p have the nice property to be 
always more efficient than M p whatever the actual statistics 
of weak lensing observables. Another issue connected to the 
efficiency of weak lensing estimators is to select the optimal 
shape for the filter Um^ ("&) which defines the aperture mass. 
Such a study is performed in Zhang et al. (2003). We shall 
not investigate this point further in this paper. 
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Table 1. Survey parameters used for the numerical results, from 
the SNAP mission as given in Refregier et al. (2004). 



Survey 


A (deg 2 ) 


rig (arcmin 2 ) 






Deep 


15 


260 


0.36 


1.31 


Wide 


300 


100 


0.31 


1.13 


Wide+ 


600 


68 


0.30 


1.07 


Wide- 


150 


150 


0.33 


1.20 



3 NUMERICAL RESULTS 
Cosmological model 

We now compute the statistics of weak lensing observables as 
they could be measured from the SNAP mission. We focus 
on the low-order cumulants of the aperture-mass and the 
shear as well as their pdf. We also consider the dispersion 
of the measures due to the intrinsic ellipticities of galaxies 
and to the finite number of cells on the sky provided by the 
survey. 

We consider a fiducial LCDM model with f2 m = 0.3, 
JIa = 0.7, Ho = 70 km/s/Mpc and as = 0.88, in agreement 
with recent observations. We shall also study in section 3.1 
the effect of small variations of these parameters onto weak- 
lensing observables. 

The many-body correlations of the matter density field 
are obtained from the minimal tree-model for the aperture- 
mass (Bernardeau & Valageas 2000; Munshi et al.2004) and 
the stellar model for the shear (Valageas et al.2004; Mun- 
shi et al.2004), coupled to the fit to the non- linear power- 
spectrum P(k) of the dark matter density fluctuations given 
by Peacock & Dodds (1996). These models are identical for 
the smoothed density field and up to the third-order mo- 
ment for any observable. From these density correlations 
one obtains all cumulants of any weak-lensing observable 
as well as its full pdf, as recalled in section 2.1. The predic- 
tions obtained in this manner have been compared in details 
with results from numerical simulations in previous works 
(Bernardeau & Valageas 2000; Valageas et al.2004; Barber 
et al.2004; Munshi et al.2004) and have been seen to provide 
good results. 



Survey properties 

Hereafter, we use the characteristics of the SNAP mission as 
given in Refregier et al.(2004), for several surveys. We recall 
their values in Table 1. The redshift distribution of galaxies 
is: 

n(z 8 ) oc z 2 e ~ (Zs/zo) and z max = 3. (65) 

The shear variance due to intrinsic ellipticities and measure- 
ment errors is <t„ = (je» j 2 } 1//2 . The survey covers an area A 
and the surface density of usable galaxies is n g . Therefore, 
we take for the number of galaxies within a circular field 
of radius 8 S : 

N = n 9 n8 2 s ^3u( U \ \ (^-Y , (66) 
V lOOarcmm /V larcmm / 

and for the number A?" c of cells of radius 8 S : 



i - 




j I iLA i i i I i i i I i i i I i i i 

2 4 6 8 10 

8, 



Figure 1. The ratio p of cosmic variance over the noise due to the 
intrinsic ellipticity dispersion, from eq. (23), for the aperture- mass 
M^p. We show the results obtained for different surveys: "Wide" 
(solid line), "Wide+" (lower dashed line), "Wide-" (upper dot- 
dashed line), "Wide>" (upper dashed line) and "Wide<" (lower 
dot-dashed line). 

Nc = = 2.7 x 10 5 ( -^—] ( — 6 -^Y 2 . (67) 

(2<9 S ) 2 v 300de g / U arcmin ,> v ; 

For the shear this number somewhat overestimates N c be- 
cause of the sensitivity of 7; s to long wavelengths, which 
would require the centres of different cells to be separated 
by more than 28 s in order to be uncorrelated. 

The SNAP mission will provide two surveys: a wide sur- 
vey designed for weak lensing with an area A = 300 deg 2 and 
a deep survey with A = 15 deg designed for the search for 
Type la supernovae. We shall refer to them as the "Wide" 
and "Deep" surveys. Following Refregier et al. (2004), in or- 
der to compare different survey strategies we shall also study 
in section 3.2 two hypothetical surveys, labeled "Wide+" 
and "Wide-" in Table 1, with the same observing time as 
the "Wide" survey (5 months) and a survey area A which 
is doubled or halved (implying a smaller or larger depth 
at fixed observing time). Finally, we also consider the two 
subsamples which can be obtained from the "Wide" sur- 
vey by dividing galaxies into two redshift bins: z 3 > z* 
(which we refer to as "Wide>") and z < 2, ("Wide<"). 
We choose z* = 1.23, which corresponds roughly to the sep- 
aration provided by the SNAP niters and which splits the 
"Wide" SNAP survey into two samples with the same num- 
ber of galaxies (hence n 3 = 50 arcmin -2 ). 

We plot in Fig. 1 the ratio p which measures the con- 
tribution to the noise from the "cosmic variance" relative 
to the effect associated with the intrinsic ellipticity disper- 
sion, from eq.(23), for the aperture-mass M ap . Of course, we 
can check that p increases with the radius S of the filter 
(i.e. with the number N of galaxies within the circular field 
of radius 8 S ). Thus, beyond a few arc-minutes the noise is 
dominated by the "cosmic variance". For the same reason, 
as compared with the fiducial "Wide" survey (solid line), p is 
larger for the "Wide-" survey (upper dot-dashed line) , which 
is narrower but deeper, and smaller for the "Wide+" survey 
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Figure 2. The variance (Mf p ) of the aperture-mass (solid curve). 
The largest error bars show the 1 — a dispersion cr(M2) from 
eq.(41) and eq.(44). The smaller error bars which are slightly 
shifted to the left show the dispersion obtained by neglecting 
non-Gaussian contributions (i.e. (M* p ) c = in eq.(41)). We also 
show the effect of a 10% increase of O m (from O m = 0.3 up 
to O m = 0.33, central dotted curve), of a 10% increase of the 
normalization erg of the density power-spectrum (from ag = 0.88 
up to ag = 0.97, upper dot-dashed curve), and of a 10% decrease 
of the characteristic redshift zo (eq. (65)) of the survey (from zo = 
1.13 down to zo = 1.02, lower dashed curve). 

(lower dashed line). On the other hand, we can see that p is 
much smaller for the low-z subsample "Wide<" (lower dot- 
dashed line), which contains fewer galaxies and has a lower 
variance (M a 2 p ), while it is almost unchanged for the high- 
z subsample "Wide>" (upper dashed line) which has also 
twice fewer galaxies but a larger variance (Mf p ). As seen in 
section 3.3, it will imply that although the skewness of the 
aperture- mass is smaller for the high- 2 subsample (and more 
difficult to measure) the deviations of the pdf V(M) from 
the Gaussian are easier to measure. The ratio p obtained for 
the shear components shows similar behaviours. 

3.1 A specific case study: Wide SNAP survey 

We first consider the case of the Wide SNAP survey, the 
properties of which are given in Table 1. We shall investigate 
the sensitivity to the survey parameters in next sections. 

3.1.1 Variance 

We show in Figs. 2-3 the variance of the aperture-mass 
(-^ap) an d of the smoothed shear component (7^) for the 
wide SNAP survey. As recalled above, this second-order mo- 
ment is smaller for the aperture-mass M ap which removes 
the contribution of long-wavelength density fluctuations to 
weak-lensing. Moreover, it bends down for small angular 
scales 9 S 1' since in this regime the projected density k(i?) 
shows more power at larger scales (i.e. the non-linear density 
power-spectrum P(k) grows more slowly than k~ 2 at high 
fc). We also display the 1 — a dispersion a(M.2) from eq.(41) 
and eq.(44) (largest error bars in the figures). The smaller 




3 4 6 8 10 



Figure 3. The variance {7^) of the smoothed shear component 
7i s (solid curve). The largest error bars show the 1 — a dispersion 
cr(M2) from eq.(47) and eq.(46). The smaller error bars which 
are slightly shifted to the left show the dispersion obtained by 
neglecting non-Gaussian contributions (i.e. ( r yf ) e = in eq.(47)). 
We also show the effect of a 10% increase of f2 m (from f2 m = 0.3 
up to n m = 0.33, central dotted curve), of a 10% increase of the 
normalization ag of the density power-spectrum (from erg = 0.88 
up to ag = 0.97, upper dot-dashed curve), and of a 10% decrease 
of the characteristic redshift zo (eq. (65)) of the survey (from zo = 
1.13 down to zo = 1.02, lower dashed curve). 



error bars which are slightly shifted to the left show the dis- 
persion obtained by neglecting non-Gaussian contributions 
(i.e. {M^ p ) c = in eq.(41)). Thus, we can see that neglecting 
non-Gaussianities slightly underestimates the dispersion of 
the measurements but the difference with the full calculation 
is rather small. The relative size of the error bars is some- 
what larger for the aperture-mass than for the shear because 
of the larger value of p. However, in both cases we can see 
that the wide survey of the SNAP mission should measure 
the variance of these weak-lensing observables up to a very 
good accuracy. We must point out, though, that this study 
does not take into account the possible systematic effects 
which might dominate the inaccuracy of the measures. 

We also display the results obtained with a 10% in- 
crease of fi m (from f2 m = 0.3 up to fi m = 0.33, central dot- 
ted curve), or a 10% increase of the normalization ag of the 
density power-spectrum (from ag — 0.88 up to ag = 0.97, 
upper dot-dashed curve), or a 10% decrease of the charac- 
teristic redshift zo (eq.(65)) of the survey (from zo = 1.13 
down to zo = 1.02, lower dashed curve). As is well-known 
and can be checked in Figs. 2-3, the amplitude of gravita- 
tional lensing distortions increases with f2 m and the matter 
density (see eq.(l)), with the amplitude ag of the density 
fluctuations (see eq.(l)) and with the redshift zo as the line- 
of-sight is more extended. As seen in the figures, there is 
a clear degeneracy between these parameters. On the other 
hand, assuming other parameters are known we see that fi m 
can be determined down to a few percents or to the relative 
accuracy of the redshift distribution and half the relative 
accuracy of ag. 

Here we must note that the different points in Figs. 2-3 
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Figure 4. The skewness 5 3 ap = (M| p >/(Mf p ) 2 of the aperture- 
mass (solid curve). The central error bars show the 1 — cr disper- 
sion (J (S3) from eq.(45). The smaller error bars which are slightly 
shifted to the left show the dispersion obtained by neglecting 
non-Gaussian contributions (i.e. (Mf p ) c = (Af* p ) c = (Aff p ) c = 
in eq. (42)). The smaller error bars which are slightly shifted to 
the right show the dispersion obtained from the estimator in 
eq. (57). We also show the effect of a 10% increase of f2 m (lower 
dotted curve), of a 10% increase of ag (central dot-dashed curve) 
and of a 10% decrease of the redshift zq (upper dashed curve). 



Figure 5. The kurtosis Sj" = (7*3)0/(7?) 3 of the smoothed 
shear component. The central error bars show the 1 — a disper- 
sion cr (1S4) from eq. (50). The smaller error bars which are slightly 
shifted to the left show the dispersion obtained by neglecting 
non-Gaussian contributions (i.e. <7f s ) c = (7f a )c = (7* s )e = 
in eq. (48)). The smaller error bars which are slightly shifted to 
the right show the dispersion obtained from the estimator S¥ in 
eq.(61). We also show the effect of a 10% increase of f2 m (lower 
dotted curve), of a 10% increase of ag (central dot-dashed curve) 
and of a 10% decrease of the redshift zq (upper dashed curve). 



are not fully independent since different wavelengths of the 
underlying density field are somewhat correlated. In order 
to combine various angular scales to obtain an overall error- 
bar on a few cosmological parameters it is convenient to 
adopt a Fisher matrix approach (e.g., Hu & Tegmark 1999). 
However, we shall not investigate this traditional approach 
here, as it has already been studied in the literature. 



3.1.2 Non-Gaussianities 

Next, we display in Figs. 4-5 the skewness S^ Ielf of the 
aperture- mass and the kurtosis Sj is of the shear component. 
These quantities provide a measure of the departure from 
Gaussianity of the underlying matter density field. They can 
also be used to break the degeneracy between the normal- 
ization of the density power-spectrum and the cosmological 
parameters (Bernardeau et al.1997). We can check that the 
error bars increase very fast with the order of the statistics. 
In particular, it is clear that the dispersion is dominated by 
the error bars associated with higher-order moments and we 
can neglect the dispersion due to the denominators (Aff p ) 2 
and (7^} 3 which enter the definition of the skewness and 
kurtosis. The central error bars in the figures show the 1 — cr 
dispersion a(S p ) from eqs.(45),(50). It is much larger for the 
shear kurtosis, which is a higher-order statistics, than for 
the aperture-mass skewness. In particular, while the detec- 
tion of non-Gaussianities from the aperture-mass should be 
clear up to 10' and even somewhat beyond, it should become 
rather difficult from the shear component for angular scales 
63 <; 6'. The smaller error bars which are slightly shifted 
to the left show the dispersion obtained by neglecting non- 
Gaussian contributions. Thus, neglecting non-Gaussianities 



again leads to an underestimate of the dispersion of the mea- 
sures, but the effect remains small. On the other hand, the 
smaller error bars which are slightly shifted to the right 
show the dispersion obtained from the estimators S p in 
eqs.(57),(61). As seen in section 2.4, these estimators which 
directly measure the cumulants always give a smaller disper- 
sion than the estimators M p which measure the moments. 
We can see in Figs. 4-5 that the improvement is rather small 
for the skewness of the aperture-mass but it is already signif- 
icant for the kurtosis of the shear. Therefore, these estima- 
tors should prove useful to extract quantitative informations 
from future weak-lensing surveys. 

As for the variance, we also display the results obtained 
with a 10% increase of f2 m (lower dotted curve), or a 10% 
increase of ag (central dot-dashed curve), or a 10% decrease 
of the redshift zq (upper dashed curve). The skewness and 
the kurtosis show a modest dependence on the cosmolog- 
ical parameter (roughly of the same order: 8%) and they 
decrease for larger f2 m . This can be understood from the 
factor fi m which appears in eq.(l). They show a somewhat 
stronger dependence on ag and they actually increase with 
ag (so that il m and ag have opposite effects on the skew- 
ness and the kurtosis while they acted in the same direction 
for the variance). This reflects the fact that a higher ag im- 
plies a matter density field which is deeper in the non-linear 
regime and further away from the Gaussian. On the other 
hand, the skewness and the kurtosis show a strong variation 
with the redshift (~ 17% and ~ 40%) and they increase 
for a smaller source redshift. This can be understood from 
the fact that a larger source redshift means a longer line-of- 
sight (whence the pdf becomes closer to a Gaussian as we 
add the lensing contributions from successive mass sheets, 



© 0000 RAS, MNRAS 000, 000-000 



10 P. Valageas et al. 




Figure 6. The pdf V{M) from eq.(26) for the estimator M as- 
sociated with the aperture-mass M ap as defined in eq. (20). Note 
that the Gaussian noise introduced by intrinsic ellipticities makes 
V(M) closer to the Gaussian than the actual pdf / P(M ap ) which 
only takes into account gravitational lensing. The solid line shows 
the theoretical prediction (26), the dashed line is the Gaussian 
and the dotted line is the Edgeworth expansion (5f) up to the 
first non-Gaussian term (the skewness). The error bars show the 
1 — <t dispersion cr(Pk) from eq.(31). 

following the central limit theorem) and a matter density 
field which is closer to Gaussian. Note that because the sen- 
sitivity onto zo and f2 m is different between the second and 
higher-order moments, they can be used to constrain both 
quantities and to remove the degeneracy between zo and 
fi m which appeared in the variance. On the other hand, as- 
suming that the redshift distribution is well known, higher 
order moments can also be used to break the degeneracy be- 
tween the normalization as of the power-spectrum and the 
cosmological parameter 57 m , as seen from the figures (also 
Bernardeau et al.1997). However, our results show that the 
error bar on the measure of f2 m cannot be smaller than twice 
the error bar on the redshift distribution. 

3.1.3 Probability distribution functions 

3.1.3.1 Aperture-mass M ap Finally, we show in Figs. 6- 
7 the pdf V(M) obtained for the estimator M defined in 
eq.(20). As seen in section 2.2, these pdfs are actually the 
convolution of the pdf "P(M ap ) (which measures gravita- 
tional lensing) by the Gaussian of variance a 2 associated 
with the noise due to the galaxy intrinsic ellipticity dis- 
persion and to the detector white noise. This convolution 
makes the final pdf V(M) closer to Gaussian than the ac- 
tual P(M B p). We consider the angular scale 9 S = 2' and 
we display the theoretical prediction (26) (solid line), the 
Gaussian (dashed line) and the Edgeworth expansion (51) 
up to the first non-Gaussian term (i.e. the skewness) (dot- 
ted line). We also show the 1 — a error bars obtained from 
eq.(31). We chose for the width A of the intervals the values 
A = (M 2 ) 1/2 /8 in Fig. 6 and A = (M 2 ) 1/2 /3 in Fig. 7. We 
can see from Fig. 6 that it should be possible to measure the 
departure from the Gaussian near the peak of V(M), which 
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Figure 7. The logarithm log('P(M)) of the pdf shown in Fig. 6. 
The error bars correspond to log(P(Af) ± a), from eq.(31). 



is slightly higher and shifted towards negative M. This also 
translates the asymmetry of V(M ap ). On the other hand, 
from Fig. 7 it appears that the far positive tail of V(M), for 
M ~ 0.05, should also provide a means to detect such non- 
Gaussianities. One should also be able to extract some useful 
information from the negative tail at M ~ —0.03. Note that 
it should be possible to distinguish V(M) from both the 
Gaussian and the Edgeworth expansion. This means that 
one has access to more information than is encoded in the 
variance and the skewness. Thus, it would be interesting to 
check in future surveys that these three domains of V(M) 
show the expected behaviour which characterizes the non- 
Gaussianities brought by non-linear gravitational clustering. 
From another point of view, the expected shape oiV{M) due 
to gravitational lensing might be useful in order to discrim- 
inate this signal from possible non-Gaussianities induced by 
the detector (which might contaminate the measure of the 
skewness) . 

3.1.3.2 Smoothed shear component 7i S We show in 
Figs. 8-9 the pdf V{Ti) obtained for the estimator Ti de- 
fined in eq.(33). As for the aperture-mass, this is actually the 
convolution of V(~/i B ) by a Gaussian of variance a 2 , which 
models the noise due to galaxy intrinsic ellipticities and de- 
tector white noise. We again display the theoretical predic- 
tion (26) (solid line), the Gaussian (dashed line) and the 
Edgeworth expansion (51) up to the first non-Gaussian term 
(i.e. the kurtosis). The 1 — a error bars are obtained from 
eq.(37) with A = (T 2 ) 1/2 /8 in Fig. 8 and A = (T 2 ) 1/2 /3 in 
Fig. 9. We see from Fig. 8 that it should again be possible 
to measure the deviations from Gaussianity near the peak 
of V(Ti) (Ti ~ 0), which is slightly higher than for a Gaus- 
sian. The non-Gaussianity might also be measured from the 
near tail of V(Ti), at Ti ~ 0.02. As for the aperture-mass, 
measuring the pdf V(Ti) in two different ranges provides 
useful information since it can be used to check the shape 
of the non-Gaussianities expected from non-linear gravita- 
tional clustering, or to discriminate the signal from spurious 
non-Gaussianities introduced by the detector. Note indeed 
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Figure 8. The pdf V(Ti) from eq.(26) for the estimator T{ as- 
sociated with the smoothed shear component 7i s as denned in 
eq. (33). Note that the Gaussian noise introduced by intrinsic el- 
lipticities makes P(Tj) closer to the Gaussian than the actual pdf 
'Pilis) which only takes into account gravitational lensing. The 
solid line shows the theoretical prediction (26), the dashed line 
is the Gaussian and the dotted line is the Edgeworth expansion 
(51) up to the first non-Gaussian term (the kurtosis). The error 
bars show the 1 — a dispersion a(Pk) from eq. (37). 
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Figure 9. The logarithm log('P(ri)) of the pdf shown in Fig. 8. 
The error bars correspond to log('P(ri) ± cr), from eq.(37). 

that one should be able to distinguish V(Ti) from both the 
Gaussian and the Edgeworth expansion. On the other hand, 
we see from Fig. 9 that the far tail of the pdf (Ft J> 0.04) 
is too noisy to give useful constraint on non-Gaussianities, 
contrary to the aperture-mass. 

3.2 Survey strategy: Width vs. Depth 

We have seen in the previous sections that the nominal 
Wide SNAP survey should yield useful information about 
the amplitude and the non-Gaussianities of the matter den- 
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Figure 10. The variance (Aff p ) of the aperture-mass for the 
Wide survey (solid line), the Deep survey (dotted line), the 
"Wide+" survey (lower dashed line) and the "Wide-" survey (up- 
per dot-dashed line). The error bars show the 1 — <r dispersion 
a{Mi) from eq.(41) and eq.(44). 
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Figure 11. The variance (7? ) of the smoothed shear component 
7i s for the Wide survey (solid line), the Deep survey (dotted line), 
the "Wide+" survey (lower dashed line) and the "Wide-" survey 
(upper dot-dashed line). The error bars show the 1 — a dispersion 
cr{M2) from eq.(41) and eq.(44). 



sity fluctuations. We now study the dependence of these 
results on the survey properties. Thus, following Refregier 
et al. (2004), we compare this Wide SNAP survey with the 
Deep survey realized by the same mission (designed for the 
search for Type la supernovae) and with two hypothetical 
surveys, labeled "Wide+" and "Wide-" in Table 1, with the 
same observing time (5 months) and a survey area A which 
is doubled or halved (implying a smaller or larger depth at 
fixed observing time). 
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Figure 12. The skewness S 3 ap of the aperture-mass for the 
Wide survey (solid line), the Deep survey (dotted line), the 
"Wide+" survey (upper dashed line) and the "Wide-" survey 
(lower dot-dashed line). The error bars show the dispersion ob- 
tained from the estimator Sg in eq.(57). 

3.2.1 Variance 

We show in Figs. 10-11 the variance of the aperture-mass 
(A^ap) and of the smoothed shear component (■yf s ) for these 
four surveys: "Wide" (solid line, as in section 3.1.1), "Deep" 
(dotted line), "Wide+" (dashed line) and "Wide-" (dot- 
dashed line). Of course, the variance is larger for the Deep 
survey since the redshift distribution is broader. However, its 
error bars are larger because the total survey area is much 
smaller. In agreement with Refregier et al. (2004), we find 
that the widest survey "Wide+" yields slightly smaller error 
bars than the nominal survey "Wide" (or the narrower sur- 
vey "Wide-" ) . Hence a wider and shallower survey is slightly 
more efficient but the difference is probably too small to have 
any impact on observational strategies. 

3.2.2 Non-Gaussianities 

Next, we display in Figs. 12-13 the skewness S 3 ap of the 
aperture-mass and the kurtosis SJ" of the shear component 
for the four surveys. The error bars correspond to the esti- 
mators H p which show less scatter than M p . The skewness 
and the kurtosis are smaller for the "Deep" and "Wide-" sur- 
veys which have a redshift distribution of sources which is 
more heavily weighted by high redshifts. As for the variance, 
the "Wide+" survey yields the best results, since it exhibits 
larger non-Gaussianities and smaller error bars. Thus, it en- 
ables one to detect non-Gaussianities up to slightly larger 
angular scales than the nominal "Wide" survey would al- 
low. 

3.3 Survey strategy: Redshift binning 

In the previous sections, we have described how the qual- 
ity of the information obtained from weak lensing measures 
depend on the survey properties. However, once a specific 
survey is realized one still has the possibility to analyze the 
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Figure 13. The kurtosis Sj 1B of the smoothed shear component 
for the Wide survey (solid line), the Deep survey (dotted line), 
the "Wide+" survey (upper dashed line) and the "Wide-" sur- 
vey (lower dot-dashed line). The error bars show the dispersion 
obtained from the estimator S^ 1 in eq.(61). 

data in different ways, for instance by subdividing the galaxy 
sample into several redshifts bins. This can be conveniently 
done by using photometric redshifts. To investigate this is- 
sue, we describe in this section the results which can be 
obtained by dividing the Wide SNAP survey, given in Ta- 
ble 1, into two redshifts bins: z s > z» (which we denote by 
"Wide>") and z < z* (which we refer to as "Wide<"). We 
choose z* = 1.23, which corresponds roughly to the sepa- 
ration provided by the SNAP filters and which splits the 
Wide SNAP survey into two samples with the same number 
of galaxies (hence n 9 = 50 arcmin -2 ). 

3. 3. 1 Variance 

We show in Figs. 14-15 the variance of the aperture-mass 
and of the smoothed shear component for the three samples: 
the full Wide SNAP survey (solid line), the high-z "Wide>" 
sample (upper dashed line) and the low-z "Wide<" sample 
(lower dot-dashed line). Of course, we find that the vari- 
ance is larger for the high-z sample since the amplitude of 
weak lensing distortions increases with the redshift of the 
source (and the length of the line-of-sight). Note that the 
error bars are quite small for all three samples therefore it 
is interesting to divide the survey into several redshift bins 
which allow one to check the evolution with time of the mat- 
ter power-spectrum. This also provides stronger constraints 
on cosmological parameters (and possibly on the equation 
of state of the dark energy, which we shall not investigate 
here) . 

We must note that the different curves in Figs. 14-15 
show some correlation since the lines of sight to distant 
sources located in different redshift bins probe the same 
density fluctuations at low z. Again it can be convenient 
to use a Fisher matrix approach to combine these vari- 
ous redshifts. We shall study the cross-correlations between 
different redshift subsamples in a future paper (Munshi & 
Valageas 2004). 
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Figure 14. The variance (Mf p ) of the aperture-mass for the 
full Wide survey (solid line), the high-z "Wide>" sample (upper 
dashed line) and the low-z "Wide<" sample (lower dot-dashed 
line). The error bars show the 1— a dispersion a(M2) from eq.(41) 
and eq.(44). 



full Wide SNAP survey (solid line), the high-z "Wide>" sample 
(lower dashed line) and the low-z "Wide<" sample (upper dot- 
dashed line). The error bars show the dispersion obtained from 
the estimator 5g in eq. (57). 
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Figure 15. The variance (7^ a ) of the smoothed shear compo- 
nent 7i a for the full Wide survey (solid line), the high-z "Wide>" 
sample (upper dashed line) and the low-z "Wide<" sample (lower 
dot-dashed line). The error bars show the 1 — a dispersion cr(A42) 
from eq. (41) and eq. (44). 



Figure 17. The kurtosis Sj" of the smoothed shear component 
for the full Wide survey (solid line), the high-z "Wide>" sample 
(lower dashed line) and the low-z "Wide<" sample (upper dot- 
dashed line). The error bars show the dispersion obtained from 
the estimator in eq.(61). 



3.3.2 Non-Gaussianities 

Next, we show in Figs. 16-17 the skewness S^ Iap of the 
aperture-mass and the kurtosis Sj is of the smoothed shear 
component for the three samples. The parameters S p are 
smaller for the high-z sample (lower dashed line) which in- 
volves the convolution of the weak lensing effects arising 
from many successive mass sheets (which makes the sig- 
nal closer to Gaussian, following the central limit theorem) 
and which probes a density field which is closer to the lin- 
ear Gaussian regime. In the case of the aperture-mass, all 
three samples allow a clear detection of non-Gaussianity and 



a rather precise measure of S^ Iap . As for the variance, it 
will be interesting to perform such a redshift binning of the 
data in order to check the evolution with redshift of S 3 ap . 
This should strengthen the constraints obtained from ob- 
servations. Moreover, it could be useful in order to discrim- 
inate the non-Gaussianities due to the non-linear gravita- 
tional dynamics from those associated with the noise which 
might be non-Gaussian (whether it comes from the galaxy 
intrinsic ellipticities or the detector itself). In the case of 
the shear component the three samples enable one to detect 
non-Gaussianity at small angular scales 8 S 4' while the 
low-z sample allows one to go up to slightly larger angles, 
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Figure 18. The logarithm log(V(M)) of the pdf V{M) from 
eq.(26) for the high-z sample "Wide>", for the estimator M as- 
sociated with the aperture-mass M ap as defined in eq.(20). The 
solid line shows the theoretical prediction (26), the dashed line 
is the Gaussian and the dotted line is the Edgeworth expansion 
(51) up to the first non-Gaussian term (the skewness). The error 
bars correspond to log('P(M) ± a), from eq.(31). 

9 s Si 7', because it yields a kurtosis which is much larger. 
Note that in both cases, the skewness or the kurtosis shows 
a strong variation with the redshift binning, which should 
easily be detected. 

3.3.3 Probability distribution functions 

Finally, we display in Fig. 18 the pdf V(M) obtained for 
the estimator M associated with the aperture-mass M ap 
for the high-z sample "Wide>". Indeed, although the skew- 
ness S 3 ap is larger for the low-z sample, its pdf V(M) is 
closer to Gaussian because the variance (M^ p ) is smaller 
so that the influence of the intrinsic galaxy ellipticities is 
larger and this turns out to be the main factor. Therefore, 
we find that for the low-z sample the pdf V(M) cannot be 
easily distinguished from the Gaussian. By contrast, as seen 
in Fig. 18 the high-z sample still allows a clear detection 
of non-Gaussianity. In fact, as for the full sample studied 
in section 3.1.3, the tails of the distribution enable one to 
distinguish V(M) from both the Gaussian and the Edge- 
worth expansion. This should again prove useful. On the 
other hand, the pdf V(Ti) cannot be distinguished from the 
Gaussian, except near its peak for the full sample. 



4 DISCUSSION 

Weak lensing surveys are already being used to constrain al- 
lowed regions of cosmological parameter space. Future sur- 
veys such as SNAP will provide a better opportunity by 
covering a large fraction of the sky. While there has been 
a tremendous progress in understanding the effect of cos- 
mological parameters on weak lensing statistics, a complete 
analysis of realistic noise contribution for various survey 
strategies is still lacking. 



In this paper we have mainly focused on realistic surveys 
such as SNAP to compute the estimator induced variance 
due to contributions from the finite catalogue size and the 
intrinsic ellipticity distribution of galaxies. Although Pois- 
son effects (due to the discrete distribution of galaxies) are 
quite small for such surveys because of the high surface den- 
sity of galaxies n g , the other contributions can play a dom- 
inant role. We study both the aperture-mass M ap and the 
smoothed shear components ji S , which can be more easily re- 
covered from actual surveys than the smoothed convergence 
k s . Extending earlier works which focused on the lower order 
cumulants, we present a unified approach in order to handle 
both cumulants and the full pdfs of these objects. 

In agreement with previous works (Refregier et al. 
2004), we find that surveys like the SNAP mission can mea- 
sure the variance of both M ap and 7is up to a very good ac- 
curacy (a few percent) for the entire range of angular scales 
1' < 9 S < 10' that we have studied. This should yield strong 
constraints on the cosmological parameters (e.g., a measure 
of f2 m to a few percent if all other parameters are known). 
However, there is a well-known degeneracy between f2 m and 
as- In addition, we find that f2 m cannot be measured to bet- 
ter than the relative accuracy of the redshift distribution, 
which might be a significant limitation. 

As usual, the degeneracy between several cosmological 
parameters can be removed by measuring higher-order mo- 
ments. We find that the skewness S3 ap of the aperture-mass 
should be easily detected and measured up to a 10% accu- 
racy over 1' < 9 S < 10'. By contrast, the shear kurtosis Sj is 
should be difficult to measure beyond 6'. Indeed, higher- 
order cumulants are increasingly difficult to measure from 
noisy data: their scatter grows with their order as a larger 
number of terms contributes to their dispersion which also 
involves all cumulants up to twice their order. Using a re- 
alistic model for the underlying matter density field, our 
computation of these error bars takes into account all these 
cumulants (that is we do not keep only the Gaussian terms or 
multiply this contribution by a fudge factor) , which slightly 
increase the dispersion. On the other hand, we introduced a 
new class of estimators 7i p designed to measure these low- 
order cumulants. We have shown that they yield a scatter 
which is always smaller than the one obtained by using the 
simple estimators M p derived from the moments themselves 
(and are actually optimal among a one-parameter family of 
estimators for a Gaussian distribution). We find that the 
gain is rather small for the skewness S 3 ap but for the kur- 
tosis SJ* S we get a significant reduction of the scatter. Since 
these estimators are no more difficult to use than the naive 
M p estimators they should be preferred over the latter. The 
skewness or the kurtosis may be used to remove the degener- 
acy between f2 m and as so as to enable one to measure the 
cosmological parameters. However, we note that they are 
very sensitive to the redshift of the sources, so that the ac- 
curacy of f2 m cannot be smaller than twice the error bar on 
the redshift of the galaxy distribution. On the other hand, 
by binning the sample over redshift one might be able to 
discriminate the influence of the redshift. 

In addition to the low-order moments we have also stud- 
ied the full pdfs V(M) and V(T X ), where M (resp. T») is 
a biased estimator of the aperture-mass (resp. of the shear 
component). Here the intrinsic galaxy ellipticity distribution 
plays a key role as it makes M and F; biased estimators and 
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it makes the pdfs V(M) and V(Yi) much closer to the Gaus- 
sian than the pdfs V{M a , p ) and V{^i B )- Note that the intrin- 
sic galaxy ellipticity also contributed to the measure of the 
low-order cumulants but it was less of a problem because 
one can still build unbiased estimators of these low-order 
cumulants (by counting each galaxy only once, which also 
removes some contributions to their scatter). We find that 
the pdfs V(M) and V(Ti) can be distinguished from a Gaus- 
sian through their shape near their maximum. Moreover, the 
negative and positive tails of the pdf V(M) associated with 
the aperture-mass can be discriminated from both the Gaus- 
sian and the Edgeworth expansion (using only the first non- 
Gaussian term defined by the skewness). This means that 
one can extract more information than is encoded in the 
first two low-order moments. Moreover, by measuring the 
pdf P(M) over these three domains one should be able to 
strengthen the constraints derived for the underlying matter 
density field and to discriminate possible non-Gaussianities 
induced by the detector. On the other hand, we find that 
the far tails of the symmetric pdf V(Ti) associated with the 
shear component are too noisy to give useful constraints. A 
detailed x 2 analysis will be presented elsewhere for simu- 
lated observations. 

Next, we have investigated whether the information ob- 
tained from observations could be improved by changing the 
survey strategy. Thus, we have compared the nominal wide 
SNAP survey with the "Deep" SNAP survey as well as with 
two hypothetical surveys with the same observing time as 
a the original survey but a different trade-off between area 
and depth. We focused on the first two low-order moments 
for the aperture- mass and the shear component. In agree- 
ment with Refregier et al. (2004) (who studied the lensing 
power-spectrum and the convergence skewness) we find that 
a wider and shallower survey is slightly more efficient. 

Finally, we have studied the possibility to divide the 
wide SNAP survey into two redshift bins (z < 1.23 and 
z > 1.23). All three samples provide a very accurate measure 
of the variance of both the aperture- mass and the shear com- 
ponent. As noticed above, this should allow one the improve 
the constraints and to check the evolution with redshift of 
non-linear gravitational clustering. The skewness S 3 ap is 
also well measured in the three samples, while the kurtosis 
Sj is is more easily obtained from the low-2 sample. This 
shows again the interest of such a redshift binning of the 
data. In addition, the high-z sample again allows a good 
measure of the pdf V(M) associated with the aperture- mass, 
which can be distinguished from both the Gaussian and the 
Edgeworth expansion. 

In this article we have mainly focused on errors asso- 
ciated with quantities derived from a specific redshift bin. 
We shall extend our studies to the cross correlations among 
various redshift bins in future works (Munshi & Valageas 
2004). This can also be generalized to compute the cross 
correlations among various statistics derived from different 
surveys with non-identical scan strategies. 

Throughout our studies we have ignored source cluster- 
ing and the effect due to lens coupling. A detailed prediction 
of source clustering and lens coupling will require an accu- 
rate picture of how galaxy number densities are related to 
the underlying mass distribution. Some of these issues have 
been studied by Bernardeau (1997), Bernardeau et al. (1997) 
and Schneider et al. (1998) who found that such corrections 



are negligible at least in the quasi- linear regime. In the non- 
linear regime one might use numerical simulations in order 
to evaluate this affect, but this would again require a specific 
recipe for the correlation between galaxies and dark matter. 

Another ingredient in our calculations has been the so 
called Born approximation. Its validity can only be checked 
numerically in the highly non-linear regime. Thus, the con- 
sistency of analytical results and numerous numerical sim- 
ulations found in previous studies strongly suggests that 
this approximation remains accurate in the highly non-linear 
regime (see also Vale & White 2003). 

In our studies we also assumed that the intrinsic el- 
lipticities of different galaxies are uncorrelated (Heymans & 
Heavens 2003, Crittenden et al. 2001) (but of course we take 
into account its variance) . Techniques to deal with such cor- 
relations have been studied in detail although the extent to 
which such correlations will affect weak lensing surveys re- 
mains somewhat uncertain. It is however generally believed 
that such effects will play a less important role as we increase 
the survey depth and we can reduce their role through acqui- 
sition of photometric redshift (Heymans & Heavens 2003). 

Although almost all present studies assume these intrin- 
sic ellipticities to be Gaussian random variables this might 
not be the case. Any signature of such non-Gaussianity if 
found by observational teams will have to be folded into an- 
alytical calculations. This can be performed in a straightfor- 
ward way within our formalism. In this case, the measure of 
full pdfs like V(M) or V(Ti) would be of great interest in or- 
der to disantengle the signal from the non-Gaussianities due 
to galaxy ellipticities (which might contaminate low-oder 
moments, especially if there are cross-correlations). How- 
ever, if such intrinsic non-Gaussianities are too large they 
might dominate the signal and preclude an accurate measure 
of the non-Gaussianities due to the matter density field. In 
addition to the galaxy intrinsic ellipticities and to the finite 
size of the survey, another source of noise is given by the 
finite PSF effects (smear) . It would be interesting to include 
this contribution in future studies. 
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